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This  comprehensive  literature  survey  reviews  past  and  present  planetary  entry  guidance  algorithms.  The 
algorithms  are  categorized  based  on  the  vehicle  L/D  and  the  planetary  atmosphere  that  the  vehicle  is 
entering.  Each  algorithm  is  described  based  on  the  guidance  type,  the  guidance  formulation,  and  the 
inclusion  of  aerothermal  parameters.  An  analysis  was  completed  defining  the  performance  of  each 
guidance  algorithm  relative  to  one  another  within  its  specific  category.  Finally,  an  overall  assessment 
was  made  regarding  the  state-of-the-art  in  spacecraft  entry  guidance. 
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1.  Introduction 

The  entry  guidance  algorithm  for  lunar  return  of  the  manned 
Apollo  capsule  (1960s)  was  the  first  flight  tested  application  of 
guidance.  Since  the  pioneering  efforts  of  Apollo  guidance 
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engineers,  many  algorithms  have  been  developed  for  different 
missions  and  vehicles.  However,  there  are  no  publications  that 
comprehensively  discuss  the  types  of  guidance  algorithms  and 
their  formulations.  This  paper  details  past  and  present  guidance 
algorithms  in  three  categories:  (1)  high  L/D,  Earth  entry  spacecraft, 
(2)  low-  to  mid-L/D,  Earth  entry  spacecraft,  and  (3)  other  planetary 
entry  spacecraft.  In  each  category,  the  vehicle  and  entry  trajectory 
type  are  defined,  see  Fig.  1  for  an  illustration  of  different  vehicle 
types.  Entry  trajectories  include  aerocapture,  direct  entry,  lofted 
entry,  and  skip  entry.  It  should  be  noted  that  the  objective  of  an 
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Nomenclature 

y 

crossrange  (nmi) 

e 

angle  between  vr  and  F  (deg) 

D 

drag  acceleration  (ft/s2) 

<t> 

declination  (deg) 

Df 

drag  force  (lbf) 

y 

flight  path  angle  (deg) 

F 

force  vector  (lbf) 

e 

longitude  (deg) 

Ft 

total  force  orthogonal  to  velocity  vector  (lbf) 

a 

bank  angle  (deg) 

Ft 

total  force  parallel  to  velocity  vector  (lbf) 

h 

planet's  rotation  vector  (1/rad) 

g 

gravitational  acceleration  [ft/s2] 

V 

azimuth  (deg) 

h 

altitude  (ft) 

ACR V 

assured  crew  return  vehicle 

L 

lift  acceleration  [ft/s2] 

AOTV 

aeroassist  orbit  transfer  vehicle 

h 

lift  force  (lbf) 

CAV 

common  aero  vehicle 

m 

vehicle  mass  (lbm) 

COV 

calculus  of  variations 

R 

magnitude  of  the  radial  distance  vector  for  3RSP 

HEMS 

high  entry  mass  systems 

model  (ft) 

MaRV 

maneuvering  re-entry  vehicle 

R 

radial  distance  vector  (ft) 

MER 

Mars  exploration  Rover-class 

t 

time  (s) 

MOrb 

Mars  2001  orbiter 

Tp 

propulsive  force  (lbf) 

MSP 

Mars  surveyor  program 

V 

relative  velocity  (ft/s) 

SC 

sphere  cone 

x 

downrange  (nmi) 

aerocapture  trajectory  is  to  meet  a  desired  orbit,  while  a  direct/ 
lofted/skip  trajectory  is  performed  to  target  a  specific  landing  site 
on  the  planet  surface,  but  both  require  a  guidance  algorithm  to  be 
meet  these  targets. 

Typically,  guidance  algorithms  developed  for  Low  L/D  ( <  0.5) 
vehicles  fly  short  range  trajectories  with  moderate  payloads 
(  —  5000  kg).  Guidance  algorithms,  developed  for  high  L/D  (>  1.0) 
spacecraft  can  fly  a  longer  cross  range,  perform  a  runway  landing, 
and  carry  very  large  payloads  to  the  ground  (upwards  of 
10,000  kg).  Guidance  algorithms  for  Mid-L/D  (>0.5  and  <  1 .0) 
vehicles  have  been  of  recent  interest  because  larger  payloads  can 
be  guided  to  the  ground  while  maintaining  low  heat  rates  and  low 
deceleration  at  higher  altitudes.  Although  the  reachable  range  for 
entry  is  dependent  on  guidance,  it  is  more  of  a  direct  function  of 
the  spacecraft  capability  and  the  entry  trajectory  type.  Some  well 
known  examples  of  low  L/D  spacecraft  are  Apollo,  Mars  Science 
Laboratory  (MSL),  and  SpaceX's  Dragon  Capsule.  The  only  high  L/D 
spacecraft  that  has  been  extensively  flight  tested  and  used  is 
NASA's  now  retired  Space  Shuttle,  whereas  most  others  are 


experimental.  Mid  L/D  spacecraft  are  currently  being  examined 
for  future  high  mass  missions  to  Mars.  There  are  other  advance¬ 
ments  in  the  works  with  the  increase  in  the  commercialization  of 
space  and  developments  in  military  hypersonics  projects.  How¬ 
ever,  many  of  these  advanced  guidance  algorithms  have  yet  to  be 
published  in  publicly  available  journals. 

This  paper  begins  with  an  introduction  to  the  field  of  entry 
guidance  and  then  discusses  the  various  models  used  to  model 
entry  flight  dynamics.  Finally,  a  comprehensive  discussion  of  the 
state-of-the  art  in  entry  guidance  is  documented  and  discussed. 


2.  Fundamentals  of  entry  guidance  development 

The  construct  for  analysis  of  past  and  present  guidance 
algorithms  is  derived  from  the  process  of  guidance  algorithm 
development.  Guidance  development  can  be  a  relatively 
subjective  process  as  there  are  numerous  ways  to  use  the  equa¬ 
tions  of  motions  and  numerous  requirements  for  a  given  mission. 


Mid  -  Law  L/D  Spacecraft 


a  -  variable  bank  angle 
a  -  fixed  angle  of  attack 


High  L/D  Spacecraft 


a  -  variable  angle  of  attack  NASA  HL-20 


Fig.  1.  Shape  and  avail 


i  controls  for  different  vehicle  shapes. 
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Trajectory  Design 


l — 

Guidance 

Development 


Trajectory  Optimization 
-classical  optimization 
-genetic  optimization 
Bank  Profile  Scan 
Determine  constraint  boi 
identify  flight  corridors 


-IE 


ConstantBank 
Modulated  Bank 
linear  ramp,  user-defined 


Reference  Trajectory  Found 


Identify’  Target: 

-Landing  Site 
-Descent  Device 
Activation  State 
-Multiple 


Trajectory  Modeling  in  Real-Time 
to  Predict  Taiget  Acquisition: 
-Profiles  Stored  in  Tables 
-Analytical  Equations 
-Numerical  Integration  of  EOMs 
-Hybrid  Approach 


Targeting  Algorithm: 
-Parameterization 
-One-Dimensional  Bank  Search 
-Optimal  Control  Formulation 
-Analytical  Equation  Blending 
-Search  Family  of  Trajectories 
-Hybrid  Approach 


Y 

|  Apply  Dispersions 
Guided  Trajectories  Found 


Fig.  2.  Trajectory  design  and  guidance  development  chart. 


The  dis/advantages  of  analytical  and  numerical  formulations. 


Analytical  formulation  Numerical  formulation 


Advantages 

Simple  to  implement  Accurate  trajectory  solutions 

Computation  time  minimal  No  simplifying  assumptions 

Solution  guaranteed 
Disadvantages 

Low  fidelity  trajectory  solution  Convergence  is  not  assured 

Formulation  tied  to  a  specific  entry  case  Extended  convergence  times 


The  flow  chart  in  Fig.  2  describes  the  process  of  re-entry  guidance 
development. 

The  first  step  to  developing  a  guidance  algorithm  is  to  design  a 
trajectory  that  takes  the  vehicle  from  the  entry  interface  to  the 
target  state  within  pre-defined  constraints  and  mission  require¬ 
ments.  This  trajectory  design  is  the  baseline  that  allows  for 
accurate  validation  of  the  nominal  guidance  algorithm.  Guidance 
algorithms  must  also  meet  the  following  in-flight  requirements 
during  the  process  of  development: 

1.  Determine  flight  feasible  control  vectors  (control  rate/accelera¬ 
tion  constraints). 

2.  Be  highly  robust  to  dispersions  and  perturbations. 

3.  Predict  trajectory  with  high  precision  and  quickly. 

4.  Acquire  target  with  high  precision  and  quickly. 

The  last  two  criteria  highlight  two  key  components  in  guidance 
development:  trajectory  prediction  and  targeting.  Trajectory  pre¬ 
diction  requires  the  developer  to  define  the  trajectory  path  with 
just  enough  detail  to  precisely  identify  whether  the  vehicle  will 
meet  its  target.  The  accuracy  of  this  prediction  is  dependent  on  the 
fidelity  of  the  formulation  relative  to  the  entry  vehicle,  desired 
trajectory,  and  available  navigational  inputs.  Note  that  during  the 
Apollo  era,  the  guidance  formulation  could  not  necessarily  use 
high  fidelity  EOMs  due  to  the  limited  computing  power  of  the 
1950s  and  1960s.  Thus,  the  Apollo  guidance  formulation  was  a 
compact  set  of  closed-form  analytical  equations  that  tracked  to  a 
reference  trajectory  which  was  stored  as  a  table.  As  computing 
power  became  exponentially  larger,  guidance  formulations  evolved 


into  numerical  solutions  integrated  with  analytical  equations, 
optimization  methods,  search  routines,  and  neural  networks.  For 
example,  the  requirements  that  NASA’s  Orion  spacecraft  complete  a 
lunar  return  in  the  contiguous  United  States  and  have  an  increased 
payload  from  Apollo  [1  ],  necessitated  the  development  of  numerical 
guidance  algorithms.  These  algorithms  had  to  achieve  extended 
range,  which  the  Apollo  Guidance  was  not  robust  enough  to 
achieve.  Despite  these  changes,  there  are  advantages  and  disadvan¬ 
tages  to  purely  analytical  and  purely  numerical  formulations,  which 
are  summarized  in  Table  1 . 

The  primary  disadvantage  to  analytical  guidance  algorithms  is 
the  application  of  several  simplifying  assumptions.  The  assump¬ 
tions  make  the  guidance  algorithm  less  robust  to  a  wider  range  of 
dispersions.  In  addition,  the  ability  to  increase  the  entry  range  or 
change  the  vehicle  aerodynamics  is  restricted  because  the 
assumptions  make  the  formulation  invalid  for  those  changes. 
The  primary  advantage  of  analytical  guidance  is  its  simplicity 
and  fast  solution  time.  Given  the  fast  timeline  for  hypersonic  entry, 
these  features  are  desirable.  Primary  disadvantages  of  numerical 
formulations  include  the  possibility  of  no  convergence  to  a 
solution  or  extended  convergence  times  to  a  solution.  This  is  not 
acceptable  since  the  hypersonic  entry  of  a  spacecraft  is  rapid  and  a 
control  solution  must  be  determined  quickly  for  a  precision 
landing. 
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Fig.  4.  Aerodynamic  force  axis  system. 


Zr 


C.G. 


r  -  Relative  Velocity 

Fig.  5.  Relative  Velocity  Coordinate  system. 

The  objective  of  the  second  key  guidance  component,  targeting, 
is  to  calculate  the  control  necessary  to  bring  the  vehicle  into  a  path 
that  meets  the  target,  without  violating  loading  constraints  and 
control  system  constraints.  The  targeting  methodology  is  also 
subjective  to  the  developer’s  viewpoint  of  what  will  be  the 
quickest,  most  efficient  method  to  find  a  new  trajectory  solution 
in  real-time.  The  flow  chart  in  Fig.  2  identifies  trajectory  prediction 
and  targeting  separately  for  general  understanding  of  guidance 
development.  However,  prediction  and  targeting  are  not  mutually 
exclusive  and  the  interaction  of  these  components  with  one 
another  is  dependent  on  the  developers  proposed  guidance 
method. 

2.J.  Entry  flight  dynamics  and  the  equations  of  motion 

The  guidance  algorithms  surveyed  in  this  paper  were,  in  part, 
analyzed  by  the  set  of  equations  used  to  model  the  predicted 
trajectory.  The  accuracy  of  the  trajectory  prediction  and  targeting 


Fig.  6.  Relative  velocity  with  respect  to  horizontal  plane. 


algorithms  depends  on  the  fidelity  of  the  models  and  the  equa¬ 
tions  of  motion  (EOMs)  it  uses.  The  following  sections  discuss 
EOMs  of  varied  fidelity. 


2.1.1.  3DOF  rotating,  spherical  planet  (3RSP)  EOMs 

The  derivation  of  this  set  of  equations  can  be  found  in  Vinh  [2], 
Note  that  the  planet  is  assumed  to  be  spherical  and  the  coordinate 
systems  in  the  two  body  problem  correspond  to  relative  motion 
dynamics.  The  Earth  coordinate  frame  is  fixed  at  the  origin  and 
rotates  with  the  Earth  (Fig.  3),  where  the  X£a^axis  is  located  in  the 
plane  of  the  Equator  and  is  pointed  toward  the  longitude  at  the 
start  of  the  trajectory.  The  Y£CF-axis  is  perpendicular  to  the  XECF- 
axis  and  the  ZEcf- axis  completes  the  right  handed  set.  There  are 
two  vehicle  axes  whose  origin  is  at  the  vehicle's  center  of  gravity. 
First,  there  is  the  axis  system  that  is  defined  with  respect  to  the 
aerodynamic  forces.  The  y-axis  is  parallel  to  the  velocity  vector,  the 
x-axis  is  parallel  to  the  vertical  lift  component,  and  the  z-axis  is 
parallel  to  the  horizontal  lift  component  (Fig.  4). 

Second,  the  relative  velocity  (r)  reference  system  describes 
the  motion  of  the  vehicle  with  respect  to  the  free  stream.  The 
positive  x-direction  is  parallel  to  the  relative  velocity  vector,  the 
positive  y-direction  points  toward  the  right  (perpendicular  to 
the  position  vector),  and  the  z-direction  completes  the  right- 
handed  system  (Fig.  5).  If  it  is  assumed  that  there  are  no  side 
forces  then  the  xrs-axis  coincides  with  the  xr-axis  and  the  bank 
angle  is  a  rotation  purely  about  the  xrs-axis. 

The  derivation  of  the  EOMs  can  now  be  carried  out,  beginning 
with  the  position  vector  (1)  illustrated  in  Fig.  7: 

R=Rir  (1) 

The  vector  R  is  parallel  to  the  xr-axis.  The  time  derivative  is  applied 
to  Eq.  (1),  which  results  in  the  following  equation: 


' dR 
dt ' 


(2) 


Note  that  the  derivative  was  taken  with  respect  to  the  coordinate 
frame  indicated  in  the  superscript  to  the  left  of  the  derivatives.  The 
relative  velocity  vector  (3)  is  a  function  of  the  speed  of  the  vehicle, 
its  flight  path  angle  relative  to  the  horizontal  plane  (Fig.  6),  and  the 
heading  angle  within  the  horizontal  plane: 


vr=V  sin  yir+V  cos  y  sin  iy]r  +  V  cos  y  cos  y/kr  (3) 

Next,  the  time  derivative  of  the  velocity  vector  is  taken  to  get  the 
relative  acceleration  vector: 


‘dvr 

dt 


-^+2 

dt2 


'M. 

dt 


+S2  x  (£2  x  R) 


(4) 


The  acceleration  equation  (4)  includes  several  terms  on  the  right 
hand  side,  the  first  term  is  the  acceleration  of  the  vehicle,  the 
second  term  is  the  Coriolis  acceleration,  and  the  third  term  is  the 
centripetal  acceleration,  all  in  the  relative  coordinate  system.  Next, 
Eq.  (4)  is  multiplied  by  the  mass  of  the  vehicle  resulting  in  the 
total  force  equation: 


T+2mQ  x  -j-r+mi 2  x  (Cl  x  R) 


(5) 
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Note  that  the  Earth's  rotation  vector  in  the  relative  coordinate 
system  is  as  follows: 

Q  =  Q  sin  (pir+C2  cos  (pkr  (6) 

The  force  vector  (F)  consists  of  gravitational,  aerodynamic,  and 
propulsive  forces.  The  gravity  force  is  acting  along  the  x-axis  in  the 
relative  coordinate  system  and  gravitational  acceleration  (g)  is 
calculated  assuming  a  spherical  Earth.  The  aerodynamic  and 
propulsive  forces  are  acting  along  and  orthogonal  to  the  velocity 
vector.  The  magnitude  of  the  force  acting  along  the  velocity  vector  is 
Ft—  Tp  cos  e-Df  (7) 


r  =  y[i-  cos  <t+  cos  +  2£2 V  cos  <p  sin  y/ 

+QlR  cos  (p(  cos  y  cos  (p+  sin  y  cos  yr  sin  (p) J  (16) 

1  \l  sin  a  V2 

yr  =  —  +—  cos  y  sin  y/  tan  <j> ~ 2QV(  tan  y  cos  y/  cos  (p 

-  sin  ip)  +  ^  y  sin  yr  sin  <p  cos  (17) 

Eqs.  (12)-(17)  are  considered  high  fidelity. 


The  magnitude  of  the  force  acting  orthogonal  to  the  velocity  vector  is 
Fn  =  TP  sin  e+Lf  (8) 

Since  Ft  is  also  acting  in  the  direction  of  the  velocity  vector,  Eq.  (3) 
can  be  used  to  write  Ft  in  relative  coordinates: 

Ft  =  Ft  sin  yirs+Ft  cos  y  sin  y/j^  +  Ft  cos  y  cos  y/krs  (9) 

The  force  orthogonal  to  the  velocity  vector  is  expressed  in  relative 

coordinates  by  applying  transformation  matrices  to  F„  in  the 
rs-coordinate  system  (10): 


2.1.2.  3DOF  non-rotating,  spherical  earth  (3NSP)  EOMs 

This  set  of  equations  follows  from  assuming  that  the  Earth's 
rotation  is  negligible.  Thus  the  EOMs  are  as  follows: 

R  =  V  sin  y  (18) 


V  cos  y  sin  yr 
R  cos  rp 


09) 


7  cos  yr 


(20) 


V  =  -D-g  sin  y 


(21) 


Fn=Fn  cos  C7zr+Fn  sin  okr  (10) 

The  bank  angle  is  a  rotation  about  the  velocity  vector,  which  is 
measured  between  the  lift  vector  and  the  xrs-axis.  The  transformation 
matrices  correspond  to  a  rotation  yi  in  the  horizontal  plane  followed 
by  a  rotation  y  in  the  vertical  plane.  The  results  are  as  follows: 

Fn  =  Fn  cos  a  cos  yirs-(FN  cos  a  sin  y  cos  y/+FN  sin  a  sin  \p)]r% 


(22) 


1  \L  sin  a 
V[  cos  y  +TT 


(23) 


Note  that  all  the  terms  containing  Earth's  rotation  have  been 
ignored. 


—  (Fjv  cos  a  sin  y  sin  y/-FN  sin  a  cos  yr)  kre  (11) 

After  considerable  manipulation  of  Eqs.  (2)— (11 ),  a  set  of  six  ordinary 
differential  equations  are  derived: 

R  =  V  sin  y  (12) 


2.1.3.  3DOF  non-rotating,  flat  earth  (3NFP)  EOMs 

This  set  of  equations  makes  the  assumption  that  the  Earth  is 
flat  in  addition  to  the  non-rotating  Earth  assumption.  In  Vinh  [3], 
the  flat  Earth  assumption  is  applied  such  that  a  slice  of  the 
planet  is  taken  at  the  Equator,  (p=  0.  Thus,  the  EOMs  now  take 
the  form: 

h  =  V  sin  y  (24) 


V  =  -D-g  sin  y+£22R  cos  (p{ sin  y 


(25) 

(26) 
(27) 


y  =  w\L  cos  ff+cos  y ( s ) 


(28) 


(29) 


2.1.4.  2DOF  Longitudinal  (2LON)  EOMs 

The  most  obvious  difference  between  this  set  of  EOMs  and  the 
above  is  the  reduced  dimensionality.  The  lateral  (crossrange)  and 
longitudinal  (downrange)  motion  are  decoupled  in  order  to 
simplify  the  problem.  The  assumption  of  a  non-rotating,  flat  Earth 
holds  for  this  set  of  EOMs 

h  —  V  sin  /  (30) 

s  =  V  cos  y  (31) 


S.N.  D'Souza,  N.  Sarigul-Klijn  /  Progress  in  Aerospace  Sciences  68  (2014)  64-74 


69 


V  =  —D—g  sin  y  (32) 


Equations  of  motion  for  planetary  entry  ordered  from  higher  to  lower  fidelity. 


3RSP  3DOF  Rotating,  spherical  planet  High  fidelity 

3NSP  3DOF  Non-rotating,  spherical  planet 

3NFP  3DOF  Non-rotating,  flat  planet 

2L0N  2DOF  Longitudinal  Low  fidelity 


Neglecting  lateral  motion  eliminated  the  azimuth  rate  equation 
and  results  in  a  differential  equation  for  the  total  range  (s  [nmi]). 
Eqs.  (30)-(33)  define  only  the  longitudinal  motion  of  the  vehicle. 


3.  Entry  guidance  survey  categories  and  definitions 

An  understanding  of  the  fundamentals  of  entry  guidance 
provides  a  window  by  which  to  analyze  the  heritage  and  modern 
guidance  algorithms  surveyed  in  this  paper.  Each  guidance  algo¬ 
rithm  was  reviewed  with  respect  to  three  characteristics: 

1.  Guidance  type  -  guidance  algorithms  can  be  of  many  different 
types,  namely: 


High  L/D  entry  spacecraft  guidance  algorith 


Test  case:  (Vehicle)  mission  Guidance  type  Guidance  formulation 


Aerothermal  inclusion 


(X-33,  L/D>  1.0)a  Refs.  [4-6] 
Direct  entry 


(X-33,  L/D  >  1.0/'  Refs.  [7,9,10] 
Direct  entry 


(X-33,  L/D  >1.0)b  Refs.  [11] 
Direct  entry 


In-flight 
reference 
generation  and 
tracking 


Numerically  solve  the  3RSP  EOMs  to  determine  feasible  reference  Heat  rate  is  a  vehicle  constraint  accounted 
trajectories  using  drag/bank  reversal  planning  (two  parameter  for  during  trajectory  planning 
search  and  mapping  techniques).  Track  reference  using 
2nd  order  tracking  law  for  linear  error 


In-flight  controls  Numerically  solve  the  3RSP  EOMs  by  incorporating  a  trajectory-  Heat  rate  is  a  constraint  included  in  the 
search/reference  segmenting  scheme  in  tandem  with  an  equilibrium  glide  condition  reference  trajectory  generation 
tracking  (2LON)  and  secant  search  method  to  determine  the  control  vector. 

A  closed-form  approximate  receding-horizon  control  law  based  on 
linearized  time-varying  dynamics.  The  EOMs  used  to  derived  the 
control  law:  Ref  [9]  -  3RSP;  Ref.  [10]  -  3NSP 


Generates  reference  trajectory  by  solving  the  2-point  boundary  The  first  phase  of  entry  incorporates 
value  reentry  problem  with  a  multi-dimensional  root  finding  a  heat  rate  tracking  law  for  a  constant 
algorithm  and  the  3RSP  EOMs  reference  heat  rate 


Direct  entry  (X-33,  L/D  >  1 ,0)b  Reference 
Refs.  [12]  Direct  entry  tracking 

(X-33,  L/D  >  1.0)c  Ref.  [13]  Reference 

Direct  entry  tracking 

(X-33,  L/Dsl.O)  Ref.  [14,15]  In-flight 

(X-34,  L/D  >  1.0)  Ref.  [27]  trajectory 

(X-40A,  L/D  >  1 ,0)d  Ref.  [21  ]  shaping/ 
Direct  entry  reference 

tracking 


(X-33,  L/Dlt  l.Of  Ref.  [8] 
Direct  entry 


Analytical 

optimization/ 

reference 

tracking 


(NASA  Space  Shuttle,  L/D  >  1.0)d  Reference 
Refs.  [17,18]  Direct  entry  tracking 


Linear  quadratic  regulator  using  2LON  EOMs 


2LON  EOMs  applied  to  range  control  using  drag  based 
approximations  w/  damped  harmonic  oscillator  control  law 
Interrogation  of  stored  PNN  database  w/  optimization 
to  shape  the  reference  trajectory.  Uses  COV  on  2LON  EOMs 


Analytical  drag  approximations  w/  optimization, 
no  specific  EOMs  needed  for  integration 


Analytical  closed  form  solutions  for  the  2LON  EOMs  that 
are  specific  to  each  phase  of  the  entry  profile 


Heat  rate  is  a  constraint  considered 
in  off-line  trajectory  optimization 
Heat  rate  constraint  in  trajectory  design 

None  specifically  stated 


Heat  rate  constraint  used  in  optimization 
and  heat  load  is  combined  with  range  in  the 
objective  function  during  optimization 


Drag  quadratic  coefficients,  for  the 
temperature  phase,  are  derived  from 


(NASA  Space  Shuttle,  L/D  >  1.0)e  In-flight 
Ref.  [26]  Direct  entry  reference 

planning  and 
tracking 
Reference 
tracking 

(Launch)  Ref.  [22,23]  In-flight 

(MaRV,  L/D  >  1.0)  Ref.  [24]  optimization 
(CAV,  L/D  >  1 .0)  Ref.  [25]  (Unsimplified) 
Direct  entry 


Generate  reference  trajectories  by  solving  an  optimization  problem  Heat  rate  and  heat  load  are  used  as  soft 
using  3NSP  EOMs  with  drag  and  lateral  acceleration  as  constraints  during  optimal  trajectory 

intermediate  controls  planning 


Numerically  solve  the  3RSP  EOMs  where  bank  modulation  Max.  Heat  Rate  and  max.  Heat  load  are 

is  a  function  of  energy  and  range  defined  as  mission  requirements 


The  Legendre  Pseudospectral  Method  is  a  hybrid  of  indirect 
and  direct  numerical  methods,  which  converts  the  optimal  control 
problem  to  an  NLP  problem  for  fast  convergence.  The  EOMs, 
constraints,  and  objectives  are  identified  by  the  user  for 
the  given  problem 


None  specifically  stated,  Ref.  [22-24],  Heat 
rate  is  included  as  a  constraint  during 
optimization  [25] 


Reference  trajectory  consists  of  energy  as  a  function  of  the 
distance-to-target  and  several  control  parameters/gains.  The 
energy  and  altitude  control  laws  are  based  on  classical  controls 


Heat  rate  and  heat  load  are  constraints  used 
during  off-line  optimization  of  the  reference 
trajectory  and  control  parameters 


a  High  robustness. 
b  Moderate  robustness. 
c  Low  robustness. 
d  Flight  tested. 
e  Managed  heat  load. 


Stored  reference  tracking  -  stored  reference  trajectory  variables 
and  determine  controller  gains  for  tracking. 

In-flight  reference  generation  and  tracking  -  represent  reference 
trajectory  analytically  or  store  in  tables,  then  modify  in-flight  to 
acquire  target. 

In-flight  controls  search  -  predict  the  trajectory  numerically  and 
perform  one-dimensional  search  to  acquire  target. 

In-flight  optimal  control  -  apply  optimization  techniques  such 
that  an  optimal  control  solution  is  found  to  null  targeting 
errors. 

2.  Guidance  formulation  -  the  formulation  can  be  analytical  and/or 
numerical  to  determine  the  vehicle  flight  path,  range  to  target, 
and  control  solution.  These  variables  are  found  using  any  of  the 
equations  of  motion  (EOMs)  given  in  Table  2. 

3.  Aerothermal  inclusion  -  indicates  whether  the  guidance  algo¬ 
rithm  includes  monitoring  and/or  management  of  heat  rate 
and/or  heat  load. 


The  guidance  algorithms  surveyed  are  detailed  in  Tables  3-5. 

Note  the  following  features  of  these  tables: 

•  References  with  specific  symbols  indicate  robustness:  ‘t'= high 
robustness,  t’= moderate  robustness,  and  ‘§’  =  low  robustness. 
The  order  of  robustness  is  based  on  the  available  comparative 
studies  to-date. 

•  References  without  symbols  may  still  be  very  robust,  but  have 
not  been  directly  compared  to  other  algorithms. 

•  References  that  include  the  symbol  indicate  guidance  algo¬ 
rithms  that  have  included  techniques  to  manage  heat  load. 

•  Algorithms  that  have  double  asterisks  (**)  indicate  that  this 
algorithm  has  been  flight-tested  (on-board  and/or  in-flight 
simulators). 


Note  that  the  authors  of  this  publication  welcome  any  feedback 
relating  to  the  technical  content  of  this  survey. 


4.  Entry  guidance  survey  results 

4.1.  High  L/D  earth  entry  vehicles 

There  are  several  different  high  L/D  vehicles  for  which  guidance 
algorithms  have  been  developed.  In  particular,  several  guidance 
algorithms  [4-15]  have  been  developed  for  the  X-33,  a  sub-orbital 
reusable  launch  vehicle  performing  a  direct  entry.  A  comparative 
robustness  study  [16]  of  these  algorithms  was  completed  and  it 
was  found  that  the  most  robust  algorithms  are  used  in-flight 
reference  generation  and  tracking  techniques.  These  guidance 
algorithms  used  a  combination  of  numerical  and  analytical  for¬ 
mulations  derived  from  the  3RSP  EOMs  (Eqs.  (12)— (17)).  Whereas 
the  least  robust  algorithms  contained  analytical  formulations 
derived  from  the  2LON  EOMs  (Eq.  (30)-(33))  and  used  reference 
tracking  only.  The  use  of  the  2LON  EOMs  assumes  decoupled 
longitudinal  and  lateral  motion.  It  is  interesting  to  note  that  the 
baseline  guidance  algorithm  for  this  comparative  study  was 
derived  from  the  Shuttle  Entry  Guidance  (SEG)  algorithm  [17,18] 
and  was  found  to  be  the  least  robust.  The  SEG  algorithm  uses  a 
reference  tracking  technique  that  adjusts  the  reference  trajectory 
in-flight  and  is  a  purely  analytical  formulation  derived  from  the 
2LON  EOMs.  Despite  the  reduced  robustness  of  the  SEG  many 
guidance  developers  take  advantage  of  the  drag  acceleration 
approach,  including  but  not  limited  to  Refs.  [4,5,19,20],  This 
approach  describes  the  reference  trajectory  in  terms  of  the  drag 
acceleration  profile,  which  is  developed  from  a  known  nominal 
trajectory. 
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Table  5 

Other  planetary  entry  spacecraft  guidance  algorithms. 


Test  case:  (vehicle)  mission 


Guidance  type  Guidance  formulation 


Aerothermal  inclusion 


(MSLab,  L/D  <  1.0)a  Ref.  [52]  Direct  entry,  Mars  Reference 
(TBB,  L/D  <  1.0)b  Ref.  [54]  Aerocapture,  Titan  tracking 
(MSLan,  L/D  <  1.0)a  Ref.  [51 )  Direct  entry.  Mars 
(MPrm,  L/D  <  1.0)b  Ref.  [50]  Aerocapture,  Mars 
(MOrb,  L/D  <  1.0)b  Ref.  [44]  Aerocapture,  Mars 
(HEMS,  L/D  <  1.0)  Ref.  [55]  Aerocapture,  Mars  Reference 
(NEV,  L/D  <  1.0)  Ref.  [53]  Aerocapture,  Neptune  tracking 
(TBB,  L/D  <  1.0)b  Ref.  [54]  Aerocapture,  Titan 
(MPrm,  L/D  <  1.0)b  Ref.  [50]  Aerocapture,  Mars 
(MSR,  L/D  <  1.0)  Ref.  [46]  Aerocapture,  Mars 
(MSR,  L/D  <  1.0)  Ref.  [45]  Direct  entry,  Mars 
(AOTV,  L/D  >  1.0)  Ref.  [41  ]  Aeroassist,  Mars 
(MSR,  L/D  <  1.0)b  Refs.  [47]  Aerocapture,  Mars  Reference 


(MSP,  L/D  <  1 ,0)c  Refs.  [48,62]  Aerocapture/ direct  In-flight 


Analytical  predictor-corrector  using  closed 
form  solutions  to  the  2LON  EOMs,  w/  reference 
trajectory  calculated  in-flight 


None  specifically  stated. 

Refs.  [41,45,46,54,50,55],  Heat 
rate  and  heat  load  are 
monitored  parameters.  Ref. 


This  algorithm  relates  the  2DOF  EOMs  to  Keplerian  motion  EOMs  None  specifically  stated 
[61,  Chapter  1  ]  to  generate  a  closed  form  solution  that  commands 
a  specific  bank  angle 


Numerically  solve  the  3RSP  EOMs  using  a  command  vi 


(MSR,  L/D  <  1.0)c  Refs.  [49]  Aerocapture,  Mar 
(MSLab,  L/D  <  1.0)  Ref.  [63]  Direct  entry,  Mars 


controls  search  includes  roll  angle  magnitudes  and  reversal  tii 


t  None  specifically  stated 


Reference  Linearization  of  the  3NSP  EOMs  u: 

tracking  a  state-dependent  Riccati  equatior 

Reference  Range  prediction  and  control  law 

tracking  polynomial  drag  approximation  as 

In-flight  Trajectory  prediction  numerically  solves  the  3RSP  EOMs  using  None  specifically  st, 

controls  search,  fourth  order  Runge-Kutta,  except  the  planet  is  assumed  to  be  an 
drag  oblate  spheroid.  Targeting  then  iterates  on  the  control(s),  jettison 

modulation  time  (and  change  in  ballistic  coefficient),  until  the  range  error  is 
within  the  tolerance.  Note  that  this  technique  is  applied  to  an 
Earth  Aerocapture  in  [58] 


erically  sc 


rs  for  specific  Heat  lo 


targeting 
None  specifically 


trigger 

stated 


a  Flight  tested. 

c  Low  robustness. 
d  Managed  heat  load. 


Scheirman  et  al.  [21  ]  applied  their  X-33  adaptive  guidance  to 
Boeing  X-40A  and  was  tested  on  the  general  dynamics  total  in¬ 
flight  simulator.  This  method  uses  in-flight  trajectory  shaping  and 
reference  tracking  techniques  to  not  only  account  for  dispersions 
but  for  control  failures  as  well.  A  database  of  trajectories  in  the 
form  of  Polynomial  Neural  Networks  (PNN)  is  compiled  off-line. 
Then  the  reference  trajectory  is  re-shaped  in-flight  by  interrogat¬ 
ing  the  PNN  and  performing  an  optimization  to  converge  to  a 
feasible  solution.  This  formulation  uses  a  combination  of  analytical 
and  numerical  techniques. 

An  in-flight  optimization  method  has  been  developed  using 
Pseudospectral  Methods  [22-25],  This  is  a  predominantly  numer¬ 
ical  method  and  has  been  theoretically  applied  to  launch  vehicles 
and  re-entry  problems  due  in  part  to  its  fast  convergence.  How¬ 
ever,  it  is  still  unclear  whether  the  fast  convergence  is  fast  enough 
for  a  planetary  entry  problem. 

Typically,  developing  guidance  for  high  L/D  spacecraft  requires 
the  engineer  to  include  some  constraint  or  control  on  the  peak 
heat  rate.  This  is  apparent  in  almost  all  of  the  guidance  algorithms 
noted  in  Table  3.  Guidance  developers  are  rarely  concerned  with 
managing  heat  load,  which  makes  sense  given  that  precision 
landing  is  of  primary  concern.  However,  there  are  two  algorithms 
which  actually  manage  the  heat  load  within  the  guidance  algorithm. 


The  first  is  Lu  and  Hanson  [8]  an  analytical  guidance  algorithm  which 
formulates  an  optimal  control  problem  that  minimizes  range  error 
and  heat  load.  The  second  is  Mease  et  al.  [26]  which  includes  a  soft 
constraint  on  heat  load  in  the  guidance  trajectory  planner. 

4.2.  Mid-  to  low-L/D  earth  entry  vehicles 

There  have  been  many  papers  published  on  guidance  algo¬ 
rithms  developed  for  the  Crew  Exploration  Vehicle  (CEV),  includ¬ 
ing  [30-34],  There  are  two  different  studies,  Refs.  [35,36],  that 
compare  the  performance  of  numerical  and  analytical  guidance 
algorithms  performing  a  skip  trajectory.  There  is  also  a  study  [37] 
that  compares  the  performance  of  guidance  algorithms  that  use 
EOMs  of  different  fidelity  for  CEV  performing  a  direct  entry.  In  Ref.  [36] 
it  was  found  that  the  adapted  Apollo  Skip  Entry  Guidance,  which  is 
purely  analytical,  is  the  least  robust  relative  to  the  numeric  predictor- 
corrector  guidance.  The  numeric  predictor-corrector  solves  the  3RSP 
EOMS  and  performs  a  one  parameter  search  to  meet  critical  points  in 
the  trajectory.  It  should  be  noted  that  this  is  only  for  entry  trajectories 
that  require  a  slap  maneuver.  In  Ref.  [35]  the  performance  of  two 
guidance  algorithms  for  CEV  was  compared.  The  primary  difference 
between  these  algorithms  was  the  event  transition  logic,  but  both  use 
a  numeric  predictor-corrector  coupled  with  the  Apollo  Final  Phase 
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logic.  It  was  found  that  the  more  robust  guidance  algorithm  used  the 
Apollo  event  transition  logic  during  a  skip  entry.  In  Ref.  [37]  the 
performance  of  the  guidance  algorithm  that  was  derived  from  the 
2LON  EOMs  was  more  robust  than  the  guidance  formulation  that  was 
derived  from  the  3RSP  EOMs  when  performing  a  direct  entry.  The 
above  comparative  studies  indicate  that  guidance  algorithms  using 
simplified  formulations  perform  better  for  mid-  to  low-L/D  spacecraft 
flying  a  direct  entry.  Whereas,  guidance  algorithms  that  use  higher 
fidelity  formulations  perform  better  for  high  L/D  spacecraft  flying  a 
direct  entry.  Finally,  if  the  mid-  to  low-L/D  spacecraft  is  flying  a  slap/ 
lofting  trajectory,  higher  fidelity  formulations  outperform  simplified 
formulations. 

Guidance  algorithm  development  was  pioneered  by  Apollo 
Guidance  engineers,  including  Moseley  [38]  and  Graves  and 
Harpold  [39],  They  had  to  develop  a  guidance  algorithm  that  took 
up  minimal  computing  resources  while  bringing  the  Apollo 
capsule,  carrying  United  States  astronauts,  safely  to  Earth.  Many 
in  the  entry  guidance  field  still  use  the  Apollo  guidance  algorithm 
given  that  it  has  been  thoroughly  flight  tested  and  is  simple  to 
implement  on-board. 

There  was  an  additional  study,  completed  by  engineers  at  NASA 
Johnson  Space  Center,  that  compared  three  guidance  algorithms 
applied  to  the  Aeroassist  Flight  Experiment  (AFE)  vehicle  [40],  The 
Analytical  Predictor  Corrector  (APC)  [41],  the  Numeric  Predictor 
Corrector  (NPC),  and  the  Energy  Controller  (EC)  guided  the  AFE 
through  an  aerocapture  trajectory.  All  three  algorithms  exhibited 
the  same  performance  in  the  presence  of  dispersions,  with  no  one 
algorithm  being  more  robust  than  the  other. 

The  European  Space  Agency  has  developed  a  guidance  algo¬ 
rithm  for  a  low  L/D  spacecraft  called  the  Atmospheric  Re-entry 
Demonstrator  (ARD)  and  a  mid  L/D  lifting  body  called  the  Inter¬ 
mediate  Experimental  Vehicle  (IXV).  The  algorithm  is  a  numeric 
predictor  corrector  that  uses  a  family  of  predicted  trajectories  to 
determine  the  final  guided  trajectory  that  nulls  the  range  error  [28], 
However,  there  is  no  comparative  study  to  determine  the  relative 
robustness  for  this  algorithm. 

The  aspect  of  including  aerothermal  parameters  in  guidance 
algorithms  for  mid-  to  low-L/D  spacecraft  is  significantly  different 
from  high  L/D  spacecraft.  Typically,  the  heat  rate  and/or  heat  load 
constraints  are  only  considered  when  designing  the  trajectory  for  a 
given  mission.  This  is  because  the  control  authority  yields  minimal 
changes  to  the  heating  profile.  The  heat  rate  and  heat  load  are 
primarily  sensitive  to  the  entry  state.  Once  the  nominal  trajectory 
is  known,  the  heat  rate  is  simply  monitored  during  active  guidance 
of  dispersed  trajectories.  In  the  event  that  trajectories  exceed  the 
peak  heat  rate  (or  heat  load)  constraint,  modifications  to  the 
vehicle  shape  and  trajectory  are  proposed.  For  low-  to  mid-L/D 
spacecraft  guidance  heat  load  constraints  are  not  considered. 


4.3.  Other  planetary  entry  vehicles 

Different  guidance  algorithms  for  Mars  and  Neptune  are 
described  in  this  section.  The  trajectory  profiles  for  other  planetary 
entry  include  aeroassist,  aerocapture,  and  direct  entry. 

A  joint  effort  between  the  Centre  National  d'Etudes  Spatiales 
(CNES,  France)  and  NASA  Johnson  Space  Center  was  conducted  to 
compare  the  Apollo  derived  Terminal  Point  Controller  (TPC)  [44], 
the  Analytic  Predictor  Corrector  (APC)  [41,45,46],  the  Energy 
Controller  (EC)  [47],  and  the  Numerical  Predictor  Corrector  (NPC) 
[48,49]  for  the  Mars  Premier  (MPrm,  L/D<  1)  vehicle  performing 
an  aerocapture  maneuver  [50],  It  was  found  that  the  accuracy  of 
the  APC,  in  the  presence  of  atmospheric  dispersions,  was  better 
relative  to  the  performance  of  the  other  algorithms.  The  more 
robust  algorithm,  with  respect  to  the  capture  corridor  and  other 
dispersions,  was  the  TPC. 


Interestingly,  the  two  guidance  algorithms  that  have  been 
improved  upon,  since  the  above  comparison  study,  have  been 
the  TPC  and  APC.  Specifically,  guidance  developers  for  Mars  Smart 
Lander  (MSLan)  [51],  Mars  Science  Laboratory  (MSLab)  [52], 
Neptune  Ellipsled  Vehicle  (NEV)  [53],  Titan  Blunt  Body  Aeroshell 
(TBB)  [54],  and  HEMS  [55]  have  used  these  algorithms.  Guidance 
engineers  for  the  Titan  aerocapture  mission  completed  a  compar¬ 
ison  of  the  APC  and  the  TPC.  It  was  found  that  the  TPC  had  a  much 
tighter  distribution  around  the  target  orbit  as  compared  to  the 
distribution  of  cases  in  the  APC.  Both  the  TPC  and  the  APC  use  the 
2LON  EOMs  to  develop  a  set  of  closed-form  analytical  solutions  to 
describe  the  reference  trajectory.  This  leads  to  very  simple  on¬ 
board  implementation,  however,  these  algorithms  demand  sig¬ 
nificant  work  to  adapt  for  different  missions  and  vehicles. 

The  TPC  algorithm  has  been  implemented  on-board  NASA 
MSLab  mission  and  is  a  purely  analytical  algorithm.  It  is  interesting 
to  note  that  the  simplicity  of  analytical  algorithms  is  utilized  in  the 
majority  of  other  planetary  entry  guidance  algorithms.  This  does 
not  mean  that  this  type  of  guidance  is  more  robust  than  numerical 
techniques,  but  that  the  predictability  of  guidance  performance  is 
high  relative  to  the  inherent  uncertainty  of  other  planetary  atmo¬ 
spheres.  There  has  been  some  development  on  numerical  type 
guidance  algorithms  for  other  planetary  entry.  Mease  et  al.  [56,57] 
have  done  significant  work  on  adapting  their  X-33  guidance  for  a 
mid-  to  low-L/D  vehicle  performing  a  Mars  direct  entry. 

There  is  another  guidance  method  for  aerocapture  and/or  direct 
entry  of  the  Mars  atmosphere,  which  modulates  drag  instead  of  the 
classical  lift  modulation  approach.  This  algorithm  uses  a  numeric 
predictor  corrector  to  solve  the  3RSP  EOMS  and  the  control  variables 
include  the  change  in  the  ballistic  coefficient  and  the  time  of  jettison 
[58,59].  Note  that  jettison  is  what  changes  the  vehicles  ballistic 
coefficient.  Preliminary  results  of  this  algorithm  show  the  potential 
to  meet  the  same  accuracy  of  current  guidance  algorithms  while 
reducing  system  complexity.  However,  a  comprehensive  study,  to 
compare  robustness  of  this  method  with  those  discussed  in  the 
previous  paragraph,  has  not  been  completed. 

The  aspect  of  including  aerothermal  parameters  in  guidance 
algorithms  for  other  planetary  entry  vehicles  of  mid-  to  low-L/D  is 
the  same  as  for  Earth  entry.  Typically,  the  heat  rate  and/or  heat 
load  are  constraints  when  designing  the  trajectory  for  a  given 
mission.  The  one  exception  is  Jits  et  al.  [60]  who  developed  a 
numerical  guidance  algorithm  for  a  Neptune  direct  entry  which 
incorporated  heat  load  as  a  trigger  to  activate  targeting.  This 
algorithm  provides  some  insight  into  the  possibility  of  managing 
heat  loads  during  entry  guidance. 


5.  Conclusion 

This  paper  identified  the  fundamentals  of  planetary  entry 
guidance  and  detailed  the  characteristics  of  the  vast  number  of 
guidance  algorithms  that  are  in  use  today.  The  comparative  studies 
for  Earth  entry  guidance  algorithms  indicate  that  simplified  for¬ 
mulations  perform  better  for  mid-  to  low-L/D  spacecraft  flying  a 
direct  entry.  Whereas,  guidance  algorithms  that  use  higher  fidelity 
formulations  perform  better  for  high  L/D  spacecraft  flying  a  direct 
entry.  Finally,  if  the  mid-  to  low-L/D  spacecraft  is  flying  a  skip / 
lofting  trajectory,  higher  fidelity  formulations  outperform  simpli¬ 
fied  formulations.  In  addition,  for  other  planetary  entries,  the 
simplicity  of  analytical  algorithms  is  preferred  given  the  relative 
uncertainties  that  exist.  This  does  not  mean  that  this  type  of 
guidance  is  more  robust  than  numerical  techniques,  but  that  the 
predictability  of  guidance  performance  is  high  relative  to  the 
inherent  uncertainty  of  other  planetary  atmospheres. 

Finally,  there  are  a  few  observations  that  can  be  made  from  this 
survey.  First,  the  effort  required  to  implement  new  missions  for  a 
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